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Strong short ranged positional correlations in- 
volving counterions can induce a net attrac- 
tive force between negatively charged strands of 
DNA, and lead to the formation of ion pairs 
in dilute ionic solutions. But the long range 
of the Coulomb interactions impedes the devel- 
opment of a simple local picture. We address 
this general problem by mapping the properties 
of a nonuniform system with Coulomb interac- 
tions onto those of a simpler system with short 
ranged intermolecular interactions in an effective 
external field that accounts for the averaged ef- 
fects of appropriately chosen long ranged and 
slowly varying components of the Coulomb in- 
teractions. The remaining short ranged compo- 
nents combine with the other molecular core in- 
teractions and strongly afftect pair correlations in 
dense or strongly coupled systems. We show that 
pair correlation functions in the effective short 
ranged system closely resemble those in the uni- 
form primitive model of ionic solutions, and il- 
lustrate the formation of ion pairs and clusters at 
low densities. The theory accurately describes de- 
tailed features of the effective attraction between 
two equally charged walls at strong coupling and 
intermediate separations of the walls. New ana- 
lytical results for the minimal coupling strength 
needed to get any attraction and for the separa- 
tion where the attractive force is a meLximum are 
presented. 

Strong Coulomb interactions in crowded nonuniform 
environments have important experimental consequences 
in a wide variety of biophysical applications ranging from 
DNA packaging in viruses to transport in ion channels 
Q II H II- They present major challenges to theory 
and computer simulations not only because of their char- 
acteristic long range but also because they can be very 
strong at short distances. We present here a new lo- 
cal molecular field (LMF) theory that averages over 
particular long ranged and slowly varying components of 
the Coulomb interactions 6] while still maintaining an 
accurate description of the short ranged components. It 
provides a general and physically suggestive theory for 
strongly coupled Coulomb systems and reduces exactly 
to the classical Poisson-Boltzmann (PB) approximation 
for dilute weakly coupled systems. 

We consider a general starting point where a molecule 
of species i, described by a rigid body frame with center 
at r^, interacts with an external field 4>fi{Yi) that comes 



from fixed charged solutes, or walls, or particular fixed 
molecules of a mobile species, as illustrated below. The 
subscript / indicates the source of the field, which we 
treat as a special fixed species /. The interaction between 
a pair of molecules of species i and j is assumed to have 
the general form Wijirij) = w^^y (r,y ) -I- Wq^ij{r^j), where 
Tij = Vj — Ti. The Ws,ij{Yij) denote general (repulsive 
core and other) short ranged intermolecular interactions. 
There are angular coordinates expressing orientations of 
the body frames that we do not denote explicitly. The 
Wq^ij{rij) arise from Coulomb interactions between rigid 
charge distributions (/^(r— r^) in the body frame of each 
molecule, so that 

where the caret denotes a Fourier transform and we as- 
sume there is a uniform dielectric constant e everywhere. 

To generate uniformly slowly varying components Ui^^j 
of the full Wq^ij = WqQ^ij + that are well suited for 
LMF averaging, we limit the magnitude of wave vec- 
tors making significant contributions to the integration 
m Eq. □ To that end we introduce a Gaussian func- 
tion parameterized by an important length scale ti that 
provides a smooth cutoff in /c-space, and write 

An An in.„\^ An, i n.„\^ ^ , ^ 

The first term on the right has all the characteristic long 
ranged Coulomb divergences as /c — > 0, but decays very 
rapidly to zero for fccr > 2. The desired slowly varying 
components arise when only this term is used in Eq. [3 
with an appropriate choice of a. For localized charge dis- 
tributions qi (k) we expand in a Taylor series about k = 
and take the lowest order multipole moment Q. This 
simplified expression defines the u^_ij we consider and 
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Separation of ionic potential 
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FIG. 1: Separation of the 1/r potential into a slowly varying 
piece erf(r/(T)/r and the short ranged remainder erfc(r/(j)/r. 
Two relevant a values are shown; a bigger a generates a more 
slowly varying long ranged component. 

thus prescribes a cr-dependent separation of the full inter- 
molecular potentials Wij{rij) = Ws,ij{Yij) + Wq^ij(vij) ~ 
Ws,zj (r»j ) + Wqo,ij (r»j ) + (r^ ) ee uq.^j (ry ) + Ul,^j (ry ) 
into short and long ranged parts. 

In r-space Eq. O becomes 1/r — eTt{r/a)/r + 
erfc(r I a) jr. Here erf and erfc are the usual error and com- 
plementary error functions. The &i:i{r/a)/r term is the 
electrostatic potential from a normalized Gaussian charge 
distribution with width a. As shown in Fig. ^ this re- 
mains smooth and slowly varying on the scale of cr, while 
decaying as 1/r at large r. This use of a Gaussian charge 
distribution is related to the Ewald sum method, which 
considers periodic images of ion configurations with em- 
bedded screening and compensating Gaussian charge dis- 
tributions. However, our focus is on the separation of the 
potential itself and not the effects of periodic boundary 
conditions and our choice of a is usually much smaller 
than that used in Ewald sum methods, which typically 
is proportional to the width of the simulation cell @ . 

The short ranged components uo,ij(ry ) define the 
intermolecular interactions in the special short ranged 
"mimic system" . They are comprised of the short ranged 
parts of the Coulomb interactions Wgo,y = Wq,ij ~ 
and the other short ranged core interactions Wg^ij'- 

ua,ijinj) = WsAjinj) + Wgo.y (r^)- (4) 

As suggested by Fig.^ a sets the range of Wqo.ij and thus 
determines an effective Coulomb core size j^] . The exter- 
nal potential from fixed charged solutes or walls (t>fi{Yi) = 
4'o.fi{^i) + 4'i.fi{^i) can similarly be separated into short 
and long ranged parts, as illustrated below. 

It is straightforward to arrive at explicit results 
for the ui_ij(r) and (j)iji{Ti). Here we consider 
localized charge distributions with a net charge 
Qi = J drqi{r) or a net dipole moment pi = J drrqi{r). 
If both molecules carry a net charge we find 



ui,ij{r) = qiqjeri{r/a)/er. The associated Coulomb core 
component is Wqo.ij = qiqjeiic{r/a)/er. Thus the results 
of Fig.n scaled by qiqj/e, give examples of possible ui^ij 
and Wqo^ij for ionic solution models [6j. For a monopole 
and a dipole we find ui_y(r) = qi{pj ■ V)[erf(r/CT)/er], 
and for dipoles ui^ij{r) — — (pi • V)(pj • V ) [erf (r/(T)/er]. 
The latter will lead to dipolar mimic systems with short 
ranged angular dependent interactions. 

Local molecular field approximation 

LMF theory introduces renormalized effective fields 
4>Rji{ri) that induce nonuniform singlet densities in the 
mimic system (denoted by the subscript R) that are sup- 
posed to equal those induced by the (j)ji{ri) in full system 
of interest: 

PR.J^{Y,) = Pf,{Y,). (5) 

This defines a general mapping relating structure in the 
mimic and full systems. Thermodynamic properties can 
be determined by integration of these structural rela- 
tions. 

We represent the effective field (t)Rji{vi) = 4)Qji{Yi) + 
4'RiJi{''^i) *s the sum of the known short ranged part 
4'o.fi{'^i) of the external field in the full system and a 
renormalized "perturbation component" (j^Rijiivi) that 
accounts for the averaged effects of the slowly varying 
interactions ui^ij. As discussed in detail in 0, by 
considering the balance of forces in the full and mimic 
systems when Eq.|51 holds, and making some physically 
motivated approximations, we find that the (t>Riji are 
determined up to a constant by the local molecular field 
equations: 

(t)RlJ^ir^) = J drj[S{f,j) + PRjjirj)] ui,j»(rji)- (6) 

Here the prime on the integral indicates an implicit sum- 
mation over all species j and an integration over the 
angles of the body frames. Long ranged interactions 
from the fixed species / are accounted for by the S{f,j) 
term, which denotes products of 5-functions describing 
the fixed location and orientation of /. 

Note that the average over the slowly varying uijiirji) 
in Eq. |H|is weighted by pRjj{Yj), the singlet density for 
species j (in the effective field of fixed species / but with 
no explicit reference to species i). This neglect of pair 
correlations between molecules at Vj and r^ characterizes 
a mean field approximation , and in most contexts this 
would represent a major source of error. However, a gen- 
eral feature illustrated by Fig. is that as a increases 
the uiji become progressively more slowly varying at 
short distances. Thus we can ensure that all the uiji 
will be slowly varying over the length scales of relevant 
local pair correlations in the system by choosing a larger 
than some state dependent minimum value crmin- This 
permits a consistent and controlled use of the mean field 
approximation in computing the average, and we antici- 
pate very accurate results from the LMF theory for any 
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c > Cniin if we properly describe the resulting density in 
the mimic system Q. 

At strong coupling we argue that crmin should be of 
order a characteristic nearest neighbor spacing a. The 
strong short ranged parts of the Coulomb interactions on 
the scale of a and below directly affect pair correlations 
between nearest neighbor molecules. These will be 
consistently described in the mimic system if we choose 
a — (Jmin of order a so that there are essentially nearest 
neighbor interactions between the effective Coulomb 
cores and the averaged effects of the ui,ji from the 
further neighbors are slowly varying on this scale. We 
illustrate below the utility of these ideas for simple 
models of ionic solutions at strong coupling. 

Size-asymmetric primitive model 

A model of great current interest is the size-asymmetric 
primitive model (SAPM) of ionic solutions, a fluid of op- 
positely charged hard spheres of different sizes in a uni- 
form dielectric continuum. The different hard sphere di- 
ameters crudely account for the different core sizes of real 
cations and anions, and there is an interesting and not 
well-understood dependence of the critical temperature 
and critical density in this model as the size or charge 
ratio is varied 9]. We consider in particular the uni- 
form equimolar system studied using Monte Carlo (MC) 
simulations by Weis and Levesque (WL), with sym- 
metric unit charges cq = qi = —q2 and a diameter ratio 
di/d2 = 0.4. Thus Ws,ij{r) = oo for r < dij = {di + dj)/2, 
and is zero otherwise, and Wq^ij(r) = qigj/er. 

WL characterize the states by two dimensionless pa- 
rameters: a reduced density p* = {Ni + N2)d2/V and 
an effective coupling strength T* = ^5/^2 (called q*^ 
in the notation of WL). Here Ib = (^e^/e is the Bjer- 
rum length, the distance where the interaction energy be- 
tween two unit charges eo equals fcsT, and (3 = (kBT)^^. 
Their simulations indicate that the critical point occurs 
at p* = 0.195 and F* = 15.15. We report results here for 
two strong coupling states with very different properties: 
a high density subcritical liquid state with p* = 1.4 and 
r* = 16, and a low density supercritical vapor state with 
p* = 0.04 and F* = 9. 

The competition between the Coulomb interactions 
and the packing arrangements of the embedded hard 
cores in the SAPM produces elaborate local structures 
in these strong coupling states as exhibited in the pair 
correlation functions gij{r), proportional to the density 
response to an external field 4>ij (r) = Wij (r) arising from 
a fixed ion of type i at the origin. Thus in LMF the- 
ory even uniform fluid correlation functions are described 
from a nonuniform point of view. These characteristic 
features can be very accurately reproduced in the mimic 
system using the strong coupling approximation (SCA). 
The SCA replaces (t>R.ij{r) by the known strong short 
ranged component 4>o^ij{r) = uo^ij(r) of the field from 
fixed ion i. This corresponds to fixing a mimic particle 
at the origin, or equivalcntly, approximating the gij{r) in 
the uniform ionic system by the go,ij{r) in the uniform 




FIG. 2: Dimensionless potentials and pair correlation func- 
tions for the SAPM at high density and strong coupling. The 
potentials use the left vertical axis. Both the full potential 
between positive ions (light grey shading) and the mimic in- 
teractions for two different values of a (darker grey shadings) 
are shown. The various pair correlation functions use the 
right vertical axis. Results for go, 11 using two different values 
of a are shown; the differences are barely visible on the scale 
of the graph. The right inset focuses on the high first peak of 
the cation-anion correlation function. 

mimic system 0. 

In Fig. El we compare correlation functions determined 
by WL for the high density state with p* = 1.4 and 
F* = 16 to MC simulations we carried out in the uni- 
form mimic system with a "molecular sized" choice of 
a — 1.2d2 using the SCA. Simulations of the long ranged 
system required careful and costly treatment of periodic 
boundary conditions using the Ewald sum method; this 
was not needed for the short ranged mimic system. De- 
spite the very different range and magnitude of the mimic 
system interactions, all the pair correlation functions are 
strikingly similar to those of the full SAPM. These func- 
tions are very different from the profiles of the associated 
hard sphere mixture with the charges set equal to zero, 
indicating the crucial importance of including the strong 
short ranged parts of the Coulomb interactions Wqo.ij in 
defining the mimic interactions in Eq. ^ Equally good 
results are found for larger values of cr, as illustrated in 
the figure, but the good agreement fails for much smaller 
tr, indicating that a-min is about 1.2^2 for this state. 

Qualitatively different structures are seen in the low 
density vapor state with p* — 0.04, and F* = 9, as illus- 
trated in Fig. 13 The simulations of WL show that oppo- 
sitely charged ions pair together with a typical spacing 
close to the minimum permitted by the hard core diame- 
ters, along with some transient formation of longer chain- 
like structures. The correlation functions between like- 
charged pairs exhibit pronounced peaks of essentially the 
same magnitude at a separation of r = 1.4c?2 = di + d2- 
This indicates the existence of small clusters of ion pairs, 
with the same peak position and amplitude for both 
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FIG. 3: Pair correlation functions for the SAPM in the low 
density ion pairing regime, with the same conventions as in 
Fig.H 

"H h" or " — \ — " configurations at the minimum dis- 
tance permitted by a linear arrangement of the embedded 
hard cores. These peaks also illustrate how counterions 
can induce an effective attraction between like charged 
objects, as discussed in detail in the next section. Very 
good agreement between full and mimic system correla- 
tion functions is achieved with a choice of CTmin — 3.0^2, 
consistent with the larger average spacing between dilute 
ion pairs. 

The clustering of the ions has probably presented the 
most severe challenges to theories of ionic systems. It is 
particularly crucial for the study of critical phenomena 
and vapor- liquid coexistence 9] . The PB approximation 
and the frequently used hypernettcd chain (HNC) 
integral equation fail to predict ion clustering; indeed 
the HNC equation has no solution in most of the ion 
pairing regime |lOj |. In contrast, the mimic system 
as described by the simple SCA already builds in the 
most important local features of ion aggregation. This 
very good agreement strongly suggests that the LMF 
theory can accurately represent the Coulomb cores that 
contribute to local correlation functions in more realistic 
models of ionic systems. Any remaining errors can be 
attributed mainly to deficiencies in the description of 
the other short ranged core interactions, thus permitting 
the efficient development of more accurate models. 

Charged walls with point counterions 

The suspension and self-assembly of highly charged 
polyelectrolytes (macroions) in the presence of mobile 
counterions (microions) is of great interest in biological 
systems 0]. These systems usually involve charge and 
size asymmetries much greater than that of the SAPM 
and are often studied by fixing a certain macroion 
configuration and computing the microion distribution 
and resulting forces on the macroions. We discuss 
here the simplest such model system consisting of 



uniformly charged infinite hard walls with neutralizing 
point counterions (and no co-ions) in a uniform dielectric 
environment. This model is simple enough that exact 
results in certain limits are known but it still 

illustrates many fundamental issues that arise from the 
interplay between long and short ranged forces in an 
explicitly nonuniform environment. It is clear from the 
previous section that the LMF theory can deal with 
more realistic models for the walls and counterions. 

One charged wall We first consider the case of a 
single hard wall with a uniform negative charge density 
at the z = plane, where we take the zero of 
electric potential energy. Without loss of generality we 
can assume that the counterions have a (unit) charge 
Co and express the wall charge density q^, = — eg/Z^ 
in terms of the length l^, of the side of a square en- 
closing that amount of charge. The potential energy 
(j)^'^{z) of a counterion at a distance z from the wall is 
27reQz/(eZ^). The Gouy- Chapman length Ic is defined 
as the distance where this potential equals fc^T, i.e., 
la = kBTellj/{2-Kel) = 1'^/{2ttIb). Ig characterizes 
the effective strength of the attractive wall-counterion 
interaction, and most counterions will be found near 
the wall in an effective slit whose width is proportional 
to Iq. Dimensionless combinations of thermodynamic 
variables in this simple system depend only on a single 
control parameter [13 ^ = Ib/Ig — ^^/(SttZq). 

As ^ increases (e.g., by decreasing T at fixed wall 
and counterion charge) counterions are driven increas- 
ingly close to the wall by the decreasing Ig- At strong 
coupling with ^ 3> 1 ov Ig Iw, most counterions 
are next to the wall and form a ("strongly-correlated") 
two dimensional (2D) liquid layer [T3. with average 
lateral spacing a ~ a2D = Iw fixed by local neutral- 
ity. There are indeed strong lateral correlations between 
the counterions in the 2D layer: the coupling strength 
Ta,„ = lB/a2D = C^/V(27r)^/^ > 1 for large ^ As dis- 
cussed above we then expect the effective Coulomb core 
size (Tniin to be of order a2D = Iw However, because of 
these repulsive cores, particles cannot stack perpendic- 
ular to the wall and remain near the narrow slit. Thus 
the density outside the slit is very low and there are only 
weak correlations normal to the wall. 

In the opposite weak coupling limit with ^ <C 1 or 3> 
Iw counterions can take advantage of the larger effective 
volume of the slit and adopt a more diffuse 3D packing 
to reduce their repulsive interactions. Crudely assuming 
all counterions are found within Ig of the wall and us- 
ing a simple cubic lattice to estimate the characteristic 
counterion spacing in this volume we now have a ~ a^u = 
(IwIg)^^^ = Zu,/(27r^)^/^. There is weak coupling between 
the counterions, with Paao = Ib/osd = i^^/^/(27r)^/^ < 1 
and here it is natural to take cTmin — Ib — Zu)(^/27r)^/^ as 
an estimate for the effective Coulomb core size Q. The 
crossover to strong coupling with essentially 2D packing 
and (Tmin — ci2D = Iw occurs for ^ of order unity, and the 
2D packing indeed provides a larger average spacing at 
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large ^. 

Quantitative results take an especially simple form 
[T^j l if we introduce a dimensionless rescaled density 
n{z/lc) = IgIwP{z) that incorporates the anticipated (f- 
dependent) scaling of the profile with Iq. Local neutrality 
requires that dzn{z) — 1, where z = z/Iq. Lengths 
scaled by Ig will generally be indicated by a tilde. More- 
over, because of the impulsive (5-function force at a hard 
wall, there is an exact relation between the pressure and 
the contact density. This yields the well-known contact 
theorem^ which implies n(0) — 1 for the contact value of 
the rescaled density at a single charged hard wall [llj . 

Exact results | l2| for n(z) are known in the limit 
^ ^ from a rigorous weak coupling expansion, which 
gives results agreeing with the PB approximation, 
npB{z) = A different strong coupling 

expansion gives exact results as ^ — > oo, nsc{z) = e~^. 
However, attempts to connect these limits by analyzing 
higher order terms in each expansion have had only 
limited success We now show that the LMF theory 
provides a simple, accurate, and unified approach for 
general ^. 

LMF equation for one charged wall We can 

take advantage of planar symmetry and integrate 
exactly over the lateral degrees of freedom in the long 
ranged parts uiji{rji) of the counterion-countcrion and 
wall-counterion interactions in Eq. |51 The resulting 
LMF equation can be written in dimensionless form for 
5i, 22 > as 

/•oo 

4>Ri{zi) = / dz2[-5{z2) + nR{z2)]G{z2,zi)- (7) 

Here 0^1 (zi) = f34'ni{zilG) and G(z2,zi) = -|zi - 
erf(|zi - 52|/ct) - 7r-i/2^e-[(^i-"^^)/'^l' + |z2|erf(|z2|/CT) + 
7r~^/^a'e~[^2/°'l is the Green's function associated with a 
normalized planar Gaussian charge distribution centered 
at Z2, with the zero of potential energy at zi = 0. 

The —6{z2) term in Eq. [3 accounts for the long ranged 
component 0i{zi) — — G(0, zi) of the dimensionless at- 
tractive potential (j)^^{zi) = zi between a counterion at 
zi and the negatively charged wall at Z2 = 0. The re- 
maining short ranged part (f)o(zi) = z\ ~ (^\(z\) of the 
wall potential is 



4>^{zx) = zierfc(zi/(T) 



(8) 



The effective field is then 0i{(zi) — 0o(^i) + 4>Ri{zi), 
with given by Eq.Qfor zi > 0, and infinity otherwise. 

mPB approximation To solve Eq. [7\ self consis- 
tently we must accurately determine the density nR{z) 
induced by /pniz). At weak coupling, neighboring ions 
interact weakly and the density response to the effective 
field is proportional to the ideal gas Boltzmann factor 
exp[— Using this approximation in Eq. we 
have a closed equation, which we refer to as the mimic 



IN, 



;=100 



Simulation 
mPB 
exp(-z) 

1/(z+bf, b~124 
no(z) 




FIG. 4: Upper grapli: (Note the log scales) Rescaled counte- 
rion density near one planar charged wall calculated by the 
full mPB theory and by the SCA compared to computer simu- 
lation data in The limiting exponential profile provides a 
good fit only near the wall. We can show analytically 7] that 
there is a crossover to an algebraic tail of the form l/{z + b)^ 
in the mPB theory for z ~ a^™ as seen in the graph. The 
vertical dashed line in both graphs indicates the value of ct^™. 
Lower graph: (Note the linear scales) The full dimensionless 
wall potential z compared to 4>r{z) and 4>oiz) from the mPB 
theory. 

Pois son- Boltzmann (mPB) equation. Moreover, we can 
show |2J for all ^ that a self consistent solution of the 
mPB equation will exactly satisfy both neutrality and 
the contact theorem. The latter implies that the density 
response takes the simple form nji{z) = exp[— </)jj(z)] 
with our choice of zero of energy. 

Remarkably, however, the mPB approximation also 
gives accurate results at strong coupling with ^ 3> 1, 
where there is an essentially 2D arrangement of the mimic 
particles in an effective narrow slit. Correlations normal 
to the wall are very weak and the Boltzmann factor again 
can accurately describe the density response to the z- 
dependent field (^p(z), as can be verified by more formal 
arguments [islllq . 

This motivates our use of the mPB approximation 
with nfl(z) = exp[-07?(z)] for all ^ in Eq. □ The mPB 
approximation is least justified at intermediate values of 
^, and it breaks down if a is chosen much larger than 
Cmin so that there would be strong direct interactions 
between further neighbors in the mimic system. Thus 
we provide a smooth interpolation between the known 
limiting values of (Tmin in the weak and strong coupling 
regimes by choosing a — cr^"" = Cmin(^B,^^), and fix 
C = 0.60 by finding the best fit to simulations 
at a moderately strong coupling state with ^ = 40. 
In this example it is numerically more convenient to 
differentiate the resulting mPB equation and solve for 
the effective force, which vanishes far from the wall, 
and then get the effective field by integration [iTj l. 
An iterative solution is straightforward and no other 
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simulation data are required. 

Results for one charged wall Figure 01 gives re- 
sults for nji{z) at strong coupling with ^ = 100. There 
is excellent agreement between the results of the mPB 
theory and MC simulations of the long ranged 

system carried out by Moreira and Netz. The log- log 
plot emphasizes that npi{z) has two characteristic 
regions. Near the wall there is an initial exponential 
decay arising mainly from particles in the 2D layer. This 
continues until about z ~ tT,iiin/2 where the density is 
very low and there is a crossover to algebraic decay as in 
the PB solution, but with a much larger effective Ig- A 
natural physical interpretation is that the small fraction 
of counterions outside the 2D layer adopt a diffuse 
PB-like profile generated by an effective wall whose 
charge density has been greatly reduced by the charge 
of the tightly bound counterions. This idea has been 
suggested before , but LMF theory is the first unified 
theory capable of describing both limiting regions and 
the crossover region. 

The density near the wall is very accurately described 
by the even simpler SCA, where ^_r(^) is approximated 
by (j)oiz). The resulting density no{z) = exp[—(f>Q{z)] 
can be written down immediately from Eq. |S1 As shown 
in Fig. ^ both (poiz) and 4>r{z) closely resemble the 
full potential zi near the wall for zi < crmin/2. But 
no (5) cannot describe the PB-like region at large z as 
does the full mPB theory and it does not obey the 
neutrality condition. This example highlights both the 
strengths and weaknesses of the SCA. When properly 
used to describe only short ranged correlations at strong 
coupling very accurate results can be found. 

Tvifo charged walls We now briefly consider two 
equally charged hard walls forming a real slit with width 
d, with neutralizing point counterions in between. At 
strong coupling and intermediate widths, the counteri- 
ons can induce an effective attractive force between the 
walls. Such effective attraction between like charged 
objects may have important experimental consequences, 
and it has generated a great deal of theoretical interest 

Reference gives an exact expression for the di- 
mensionless (osmotic) pressure P = PIqI^P arising from 
neutralizing point counterions confined between charged 
hard walls at z = and z = d: 

P = n{0) - 1. (9) 

Thus if the rescaled contact density n(0) is less than 
(greater than) one there is an effective attractive (re- 
pulsive) force on the walls. As d oo, we recover 
the one- wall results discussed earlier, where P = 0, and 
n(0) = 1. 

Since the total force on a counterion from equally 
charged walls exactly cancels for all z and all cZ, we now 
have 4P''^{z) = for < z < d. As in the one- wall case. 



it is useful to divide into a short ranged part 

0g»(i; d) ^ 01"' (Z) + (J- ~Z) - ~^l^{d), (10) 

given by a sum of short ranged single-wall terms (indi- 
cated by the superscript Iw) defined in Eq. |S1 and the 
remainder. We take the zero of energy on the left wall at 
z = 0. The effective field 4>b{z) = 4>'o"{z; d)+4>Ri{z) is de- 
termined from the two-wall LMF equation. This closely 
resembles Eq. [7| except that the integration is from 
to d and there is an additional —S(d — Z2) term in the 
integrand, accounting for interactions with the second 
wall aX z — d. Again we use the mPB approximation 
nR{z) = Aexp[— 0jj(z)] and fix the constant A (which 
equals the contact density with our choice of zero of en- 
ergy) using the neutrality condition dznii{z) =2. The 
pressure is then given by Eq. |51 

The resulting two-wall mPB equation reduces exactly 
to an integrated form of the PB equation if cr = 0. The 
latter has an analytic solution and predicts a repulsive 
force for all d and f • The mPB theory also predicts 

a weak repulsive force at strong coupling for sufficiently 
large d, arising from weak repulsions between counteri- 
ons in the dilute PB-like tails of the one-wall profiles 
discussed above. On the other hand, at strong coupling 
and sufficiently small d, core repulsions make it unfavor- 
able for particles in the narrow slit to stack perpendicular 
to the walls, and counterions will be forced into a single 
2D layer with characteristic lateral spacing d ±2 /^/\/2 
fixed by neutrality. To interpolate between this limit and 
weak couphng we choose a — = C mm{lB,lw/V^), 
and take the same value of C = 0.60 that we used for 
the one-wall theory. Since (Tmin — at small ^, the PB 
approximation is consistent only as ^ ^ 0. The mPB 
theory naturally introduces a crucial new length scale 
CTniin that allows for a change in the functional form of 
riR^z) as ^ increases. 

Figure compares numerical results of mPB theory 
[T7I I to simulation data Il3| for strong coupling states. 
The left graph shows for ^ = 100 there is very good agree- 
ment between the mPB theory and computer simulations 
for all widths where simulations can be performed. As 
shown in the inset, the mPB theory predicts that at still 
larger widths there is a weak repulsive force between the 
walls. This reentrant behavior is illustrated more gener- 
ally in the right graph, which also shows that a minimal 
coupling strength of ^ > — 12 is needed to get any 
attraction |T^. 

The left graph also shows results from the ana- 
lytic SCA, where n{z) is approximated by no (-2) = 

exp[— ^o"" (-^i '^)]i with Ao similarly determined by neu- 
trality. At strong coupling and large d, 4)'^^{z]d) has a 
deep attractive well near each wall essentially identical to 
that found near a single wall. However at small enough 
d the wells from the individual ^J*" terms in Eq. ^] be- 
gin to overlap and their depth decreases. As d ^ 0, 
(p^iz^d) vanishes for all z in the slit. The SCA fails 
at large widths, just as it did far from the wall in the 
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FIG. 5: Left graph: Dimensionless pressure at strong cou- 
pling between two equally charged hard walls as a function of 
width d from [T^ compared to predictions from the full mPB 
theory and the SCA. The inset shows that the mPB theory 
predicts a weak repulsive force at larger widths. Right graph: 
Repulsive and attractive forces between two charged walls as 
as a function of width and coupling strength as determined 
by MC simulations and the mPB theory. 

one-wall case, and the reentrant behavior is completely 
missed. However, it is very accurate at smaller widths 
and it gives a good description of the location and mag- 
nitude of the maximum attractive force. 

The formation of the single 2D layer at sufficiently 
small widths d < (t^™ plays a key role in producing 
a strong attractive force. Because of the absence of 
correlations normal to the walls, the density profile will 
be relatively constant across the slit. We can use the 
SCA to describe several features analytically in this 
regime. The minimum pressure P* should occur near 
the largest width d* = 2z^ for which the single 2D layer 
remains stable, defined by 0o'"(zj^;d*) = 1. At larger 
widths there will be higher contact densities as separate 
2D layers at each wall begin to form, and the nearly 
constant profiles at smaller widths will have higher con- 
tact densities by normalization. Expanding (j)Q'^{z;d) in 
a Taylor series, we have (i)l'^{z;d) = 2z{d— z) / {ir^/'^a^'") 
to lowest order; higher order terms are negligible for 
J < CT^™. This implies d* = (27ri/25-2™ )i/2 ^ 1.94^i/4^ 

using our expression for (t^™. Similarly evaluating Aq, 
the minimum pressure is P* ~ 3.717/d* — 1. When 
P* = there can be no attractive forces; this determines 
the minimal coupling strength — 13.48 and associated 
critical spacing d* ~ 3.717 needed to get any attraction. 
Finally, for ^ > and d = 2 we see 4>l'"{z]d) <C 1, 
so the constant profile nQ{z) = 2/c? is very accurate. 
Equation |51 then implies that at very strong coupling 
the transition from strong repulsive to str ong attractive 
forces occurs near d = 2, as shown in [T^ Il3j| . All these 
predictions are in very good agreement with numerical 
solutions of the mPB theory (and with simulations of 



the full and mimic systems, where available) for all 
strong couphng states, as illustrated in Fig. 

Discussion 

In strong coupling regimes the short ranged parts of the 
Coulomb interactions efficiently compete with the other 
short ranged molecular core interactions and strongly in- 
fluence pair correlations between neighboring molecules. 
In LMF theory the choice of (Xniin determines the strength 
and range of these important Coulomb core interactions. 
They play a key role in inducing an effective attraction 
between like-charged objects at strong coupling, as il- 
lustrated here in Fig. O for the correlation functions be- 
tween like-charged ions in the SAPM, and in Fig. [S] for 
the osmotic pressure on two like-charged walls. In both 
cases, strong short ranged forces mainly involving single 
counterions or a single counterion layer can mediate a 
strong effective attraction. Such phenomena have tradi- 
tionally been interpreted as illustrating the "breakdown 
of mean field theory" and the need for new and more 
sophisticated approaches. But LMF theory using the 
simple strong coupling approximation, where only the 
short ranged parts of the Coulomb interactions are taken 
into account, provides very accurate results at small and 
moderate separations. At weak coupling and for long 
wavelength correlations the averaged effects of the long 
ranged interactions as determined by the LMF equation 
are needed as well 0. 

LMF theory provides a general conceptual framework 
that connects and clarifies previous work on systems with 
both short and long ranged interactions. It has been used 
to describe liquid-vapor interfaces and wetting and dry- 
ing transitions for simple fiuids and hydrophobic in- 
teractions in water It suggests new and simplified 
simulation models for general Coulomb systems based on 
the mimic system that do not require special treatment 
of periodic boundary conditions. The short ranged inter- 
molecular interactions in the Coulomb mimic system are 
reminiscent of the truncated interactions used in reaction 
field methods . But it was not very clear in those ap- 
proaches how to choose an appropriate cutoff and how 
to treat nonuniform environments. The determination 
of (Tmin and the effective field in LMF theory provides a 
way to deal with both problems. We believe the LMF 
picture will prove useful not only in formal theory, but 
also for qualitative reasoning and in detailed simulations 
of biophysical systems. 

Detailed results for charged hard walls, including MC 
simulations of the nonuniform mimic system using both 
the SCA and the full LMF theory, with important contri- 
butions made by Charanbir Kaur and Jocelyn Rodgers, 
will be presented elsewhere |T3|. We are grateful to 
them and to Michael Fisher for many helpful remarks. 
We thank J. J. Wcis for sending us his simulation data 
for the SAPM correlation functions. This work was 
supported by the National Science Foundation through 
grants CHEOl-11104 and CHE05-17818. 
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